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Abstract. Rotation curves of spiral galaxies are known with reasonable precision for a large number of galaxies with 
similar morphologies. The data implies that non-Keplerian fall-off is seen. This implies that (i) large amounts of dark 
■ matter must exist at galactic scales or (ii) that Newtonian gravity must somehow be corrected. We present a method 
, for inverting the integral relation between an elemental law of gravity (such as Newton's) and the gravitational field 
generated by a thin disk distribution with exponential density. This allows us to identify, directly from observations, 
extensions of Newtonian gravity with the property of fitting a large class of rotation curves. The modification is inferred 
from the observed rotation curve and is finally written in terms of Newton's constant or the effective potential of a 
test mass moving in the field generated by a point-like particle. 
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When trying to understand the dynamics of large scale astrophysical systems, for example galaxies, we assume 
^ that the dominant interaction at that scale is gravity. This implies the use of Einstein's General theory of Relativity 
5^ ■ (GR), so that in the limit when the speeds involved are much smaller than the speed of light and in the weak field 
"qJj \ limit, one may legitimately apply the Newtonian limit. For most galaxies these two conditions are met, and therefore 
d • Newtonian considerations apply. 

Both GR and its Newtonian limit, have been successfully and directly tested at scales not much larger than the 



Solar System (See, e.g., Will 1993). However, when one tries to apply them to galaxies or even larger systems, the 



predicted behavior is usually found to be quite different from what is observed. In fact, in order to accommodate the 
■ observations, it is customary to assume the presence of a large amount of non- visible matter, the so-called dark matter. 
The needed amount of dark matter has to be somewhere between 90 and 99 per cent of the total mass of the Universe; 
furthermore, in order to be consistent with the predictions derived from standard Big Bang nucleosynthesis, it must 
be non-baryonic. Of course this discrepancy is commonly known as the Dark Matter Problem. 

However, Newton's law of gravity is just a phenomenological law that was designed by Newton to explain gravita- 
tional dynamics within Solar-system scales. On the other hand, General Relativity (GR) was developed by Einstein 
with the constraint in mind of recovering the Newtonian potential in the limits of weak fields and small (compared 
with the speed of light) velocities, assuming that the phenomenological law discovered by Newton for the Solar system 
could be extrapolated to distances up to infinity. These considerations leave open the possibility that, at least while 
the Dark Matter component remains unidentified, the possibility exists that GR, despite its conceptual beauty, would 
have to be modified in some way, perhaps in the same spirit as it was used to modify Newtonian gravity for strong 
fields and large velocities or, perhaps, in other ways. 

In this paper we study the problem of the rotation curves of spiral galaxies, a case in point. When we apply the 
Newtonian approximation to these systems we find that their rotation curves should fall for large radius as v 2 ex r , 
i.e., in a Keplerian fall-off. Instead, the observed rotation velocity is typically seen to remain constant after attaining 
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a maximum value, as if it were to go to some asymptotic value, different for each galaxy. This is usually explained 
assuming a halo of Dark matter surrounding the visible galaxy, with the adequate shape for accommodating the 
observed rotation curve. 

We have posed the following (somewhat longish) question: "Is it possible to find a phenomenological universal 
Newton-like law that can explain the observed dynamics in spiral galaxies, without having to assume the presence of 
an undetected mass component, and whose short distance limit is compatible with the Newtonian law?' We will offer 
a positive answer to this question. 

This is not the first time such an approach is taken in the literature. Work along these lines has already been done 



(for example by Kuhn & Kruglyak 1987 and Tohline 1983), and assumes a specific functional form for the generalized 



force. This form is parametrized by some free parameters, and the evaluation of the predicted rotation curves for 
some known spiral galaxies through the corresponding numerical integrals leads, when compared with observations, to 
specific values for these parameters. This is a very interesting and direct approach, but its success obviously depends 
on a good choice for the initial of for the "improved" force. 

We will present here a procedure that follows the inverse methodology: we will write down an equation such that, 
once we know the observed rotation velocity of a galaxy we readily obtain which is the force, if any, that is able to 
generate that rotation curve without assuming the presence of any dark matter. In this way we will not have to assume 
a form for the phenomenological law, we will infer it directly from the observational data. The observed data is our 
starting point, not the final result of some "fit" . And, what could be more interesting, the method presented here can 
be, in principle, equally useful for discarding a non-Newtonian law of gravity as for proving its existence. Once we 
have the equation that allows us to find the force from the observed velocity, we will apply it to a sample of spiral 
galaxies and check if there exists a common phenomenological law that works for all the galaxies in the sample. 

We write the gravitational field as a generalization of both, the Newtonian potential and the Newtonian force. This 
we do by introducing two functions g(r) and g c s{r) defined as: 

_ G mim 2 , . * , 

<\>{r) = g(r). (1) 



G mim 2 , , 
F W = ~2 9cff{r)- (2) 

where <p(r) and F{r) are the potential and the force experienced by two point-like particles of masses mj and 777.2 
separated by a distance r and Go is Newton's constant. The two functions g(r) and g c s(r) are related through: 

g eS (r) ee g(r) - rg'(r), (3) 

where the prime denotes a derivative with respect to the argument of the function. 

In order to calculate the field due to a given mass distribution Q described by a density function p(r), we must 
first integrate the microscopic field over the distribution, i.e., we must perform the integral 

which gives the potential experienced by a point mass at a position R from the center of O. For a symmetric distribution 
(spherical or disk when considering the disk plane) it is convenient to introduce the following notation: 

m = - ^ w (5) 



V? ot {R) ee * eff(j R) (6) 

where V ro t(R) is the rotation velocity of a test particle describing circular orbits in the gravitational field generated 
by the distribution f2. 

It can be readily checked that the auxiliary functions and W satisfy the following functional relationship: 



* cf f (i?) = - m'(R). 



(7) 
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As we have pointed out above, the luminous matter in many spiral galaxies can be well modelled by a thin disk 



distribution with an exponentially decaying density (see Freeman 1970). In cylindrical coordinates this may be written 
as: 



p(r) = po e~ ar 5(z) 



(8) 



where a is obtained from the luminosity profile for each galaxy. 

Our problem can now be paraphrased as follows: "knowing the rotation velocity (i.e., ^ c ff (R) up to a normalization 
constant proportional to the mass of the galaxy), what is the microscopic law of gravity, i.e., g(r) or g e s(r), capable 
of generating that velocity field in a thin disk galaxy?" 

This problem can be solved exactly for a spherical galaxy with an exponentially decreasing density. Here the 
solution can be summarized as: 



g (r) = *( r ) 



^'(r) + i$C)(r). 



(9) 



We will consider this case in detail in a separate publication (See Rodrigo-Blanco 1996a). The line of reasoning 
leading to the proof of Eq. (|]) is as follows: plug Eq. ^ into (|4|), and use the Fourier transform of g(r) together with 
the addition theorem of Bessel functions to decouple the integration variables in the integrals. Once this is done, and 
after introducing the function ^(R), integration by parts yields Eq. (^). 

In the thin disk case the problem cannot be solved exactly. For this reason we will use an approximation that 
we call 11 Gaussian approximation" (we will see that, in the Newtonian case, this approximation is equivalent to using 
Gauss' law for calculating the gravitational field and hence the name). It should be noted that this approximation 
improves when one considers a g c s(r) which is an increasing function of r (which, of course, is a welcome bonus for 
understanding the rotation curves of galaxies) . 

In this approximation the appropiate g e s(f) turns out to be: 



Seff(z) = ^Mx) - -Wrfr(aO + 4-*eff(z) 
or orx 

where \E , e ff(a;) has the following behavior at the origin: 
* cff (0) = ^ ff (0) = 0. 



(10) 



(11) 



The mathematical formalism applied to obtain Eqs. (|10| ) and (|ll| ) is very similar to the one used for a spherical 
distribution. In both cases, the use of the addition theorem of Bessel functions leads to an infinite series of terms 
involving Bessel functions of the form J2k+i/2 an d we truncate the series keeping only the term with k = 0. Actually, 
in the presence of spherical symmetry this is the only term that contributes to the integrals, and therefore the result 
is exact. In the thin-disk case it can be shown that this term dominates the in tegrals in the cases of interest. The 
mathematical details will be given in a separate publication (See Rodrigo-Blanco 1996b ). Here we give a qualitative a 
posteriori justification of the goodness of the approximation. First, it can be seen that the solution to Eq. (|Io|), when 
<?efif(?") = lj (Newtonian limit) is 



G M U 



[1 - (1 + uR)e- aH ] 



G M(R) 



(12) 



R ' R ' 

where M (R) is the disk mass inside a sphereB of radius R. In order to get an idea of what happens for a growing g c s , 
let us restrict ourselves to the case when g e s(r) can be parametrized as a power law of the form g^(r) = (^) M with 
/i real and positive. In Fig. (Q) we have plotted the rotation curve obtained from Eq. (|l0|) versus the exact solution 
for six values of p. It can be seen right away that, the faster g^g grows, the better the approximation. 

Now we move on to apply our equation (Eq. (|l0|)) to real galaxies. In order to do this we have chosen a set of 9 
spiral galaxies whose luminosity profiles can be well fit ted using a thin disk model with e xpone ntial density as the one 
assumed to obtai n Eq. ( |Io| ) (See refs. Begeman 1988 , Carignan 1985 , Persic & Salucci 1995 , Persic et al. 1995 and 
Mathewson et al. 1992 ). In Table |l| we list the relevant observational data for these galaxies. 

The first step for applying Eq. ([To]) is to fit the rotation velocity for each galaxy by some function, so that we 
can take its derivatives. In order to do this it is important to notice that, in the approach we are describing, there is 



1 Although this is not the exact result, it is however what we would find if we applied Gauss' law as an approximation for 
evaluating the gravitational field. That is why we call our approximation Gaussian. 
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Fig. 1. Exact rotation curves obtained by performing a numerical integration of the forces over a flat disk and obtaining the 
rotation velocity (full points) and Gaussian approximation (solid line) for g c g(r) oc r M for some values of /j, (// = is the 
Newtonian case). In every case, for the sake of clarity, the velocities are normalized by dividing by the appropriate constant: 
Vb = G °^fa)£ a • ^ n ^ ne inset graphs (dashed lines) we have plotted for each case the percentage of error made when we use the 
Gaussian approximation instead of the numerical integrals, as a function of r. 
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Fig. 2. Fits of the rotation curve of NGC 6503 using ten different functional forms (We have used the class of functions vi, 
defined in Eq. (13), with a third-degree polynomial and twelve different values of fi going from 1.0 to 2.2 (upper graph) and the 
g e «{r) corresponding to each fit (lower graph). 
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no specific physical reason or prejudice for choosing one function or another to fit the observed data: we only need 
that the function fits well the observed velocities within the error bars. Apart from this generic requirement, we can 
use powers, polynomials or any other function that fits the data so long as the corresponding velocity satisfies the 
following sufficient conditions: v(r = 0) = 0, and v(r)e~ ar — > as r — > oo (which are both natural conditions). It is 
also important to notice that, since second derivatives of the rotation velocity are involved in Eq. ([l0]), small differences 
in velocities can still lead to large differences in g c g. This can be readily appreciated in Fig. |2|, where we have plotted 
different functional fits to the rotation velocity of NGC 6503 and the corresponding g e g 's (up to a normalization 
constant). 

Although all the fits are statistically acceptable, we see qualitative differences between the different functions used 
to fit and then represent the data points. The differences are more significant at short distances. This range of distances 
can be seen in Fig. |l| to be the one for which our approximation is worse. Thus we will discard any values in those 
regions. Actually since these points are in any case out of the range where our method applies, we are justified in 
removing them since we are mainly interested in the behavior of g c g(r) for large r; the only restriction is that they be 
compatible with g e g = 1 at small distances. 
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Fig. 3. Observed rotation curve vs. the rotation curve generated by the g e « (r) selected for each galaxy. 



Table 1. Relevant observed properties of selected galaxies. References to the original data: (a) Begeman 1988; (b) Carignan 
1985; (c) Carignan et al. 1988; (d) Carignan & Puche 1990; (e) Kent 1987; (f) Mathewson et al. 1992; (g) Metcalfe & Sanks 
1991; (h) Persic & Salucci 1995, Persic et al. 1995; (i) Puche et al. 1990 



Galaxy name 


Distance 


Scale length 


Luminosity 


Rot. 




(Mpc) 


(kpc- 1 ) 


10 10 L o 


curve 


NGC 2403 


3.2 a 


2.1 a 


0.8 a 


a 


NGC 3198 


9.4 a 


2.4 e 


0.9 a 


a 


NGC 0598 


0.9 9 


im h 


0.36 h 


h 


NGC 6503 


5.9 a 


I. IT 


0.48 a 


a 


NGC 0247 


2.5 b 


2.9 b 


0.24 b 


d 


NGC 0300 


1.9 6 


2.0 b 


0.24 b 


i 


347-g33 


20.9' 


1.46 h 


1.675 h 


h 


UGC 2259 


9.8 c 


1.34 e 


0.1 c 


d 


NGC 7339 


20.6 h 


1.9 h 


1.159 h 


h 
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In view of the above we use the following procedure: (i) we fit each rotation curve by a wide family of mathematical 
functions, (ii) we calculate the corresponding product g e g(r) Mi for each of those functions (where i denotes the 
particular galaxy). Thus we have a set of g e ff(r) 's for each galaxy; furthermore, each of them can generate the 
observed rotation curve up to some multiplicative constant. Then, (m) we introduce all the g c s{f) 's in a computer 
program that picks up a g c s{r) Mi for each galaxy in such a way that, once divided by an appropriate constant, all the 
ges(f) 's are as similar as possible. For doing this, we choose a standard galaxy among the ones in our sample (in this 
case we have chosen NGC 6503 because the range of distances for which its rotation curve is observed is the best one 
to compare with the other galaxies in the sample). Then, for each galaxy, we fit the mass proportionality constant for 
each g c g(r) minimizing the \ 2 of the comparison with one g c s(f) for NGC 6503. Finally we pick, for each galaxy, the 
g e s{i") for which the final \ 2 is smallest. In this way we obtain a g e s(r) for each galaxy that is capable of generating 
its observed rotation curve within the observational accuracy, and we also obtain the total mass of the galaxy relative 
to the mass of NGC 6503. 
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Fig. 4. Normalized g c g(r) for all the galaxies in the sample. All of them are multiplied by a common constant factor 
£ = M 65 o 3 /lO lo M . 

Although, as mentioned before, it is irrelevant what class of mathematical functions we use for the fit, it is 
nevertheless interesting to mention what functions we have used here. We have used two kinds of combinations 
between powers and polynomials, labelled as v± and xi2, and defined by: 



P m (r) 



(13) 



v 2 2 {r)=r t P m {r) 



(14) 



where, in each case, P m (r) is a polynomial of degree m in r and t is an integer greater than or equal to one. 

In Fig. [3] we show the fit to the rotation curve for each galaxy in our sample, and in Fig. [| we plot the corresponding 
g e s(r) for the galaxies multiplied by the mass of NGC 6503 in units of 10 10 solar masses. In Table |^ we list the mass 
of each galaxy in terms of the mass of NGC 6503 as well as the corresponding mass to light ratio for each galaxy. This 
mass-to-light ratio is in units of Mq/Z b q and Mq^/I^Mq. In this table we also indicate which type of function 
V\ or V2 (See Eq. @ and Eq. @) was finally chosen to fit the observed rotation curve of each galaxy. 

We end by offering some conclusions. 

We have found the solution to the problem of inverting the integral relation between an elemental (two-body) law 
of gravity and the gravitational field generated by a thin disk distribution with an exponentially decaying density. We 
have solved this problem in an approximation that we have called Gaussian. We have shown that this approximation 
in general leads to good results at large distances, although it fails at short distances where (in any event) Many Body 
effects may be relevant and overshadow the physics of few bodies. This, together with the facts that rotation curves 
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Table 2. The first column indicates which kind of function vi t>2 (See text) was chosen to fit the observed rotation curve for 
each galaxy. The second and third columns show the parameters \i and m that better fit the data. In the fourth column the 
mass corresponding to the g e s (r) chosen for each galaxy is given. The last column shows the corresponding mass-to-light ratio 
calculated using the observed value of Lb this mass-to-light ratio is in units of Mq/L b q and Mq$o3 / 10 10 Mq. 



Galaxy 


t 1 


m 




M/M 6503 


(M/L B ) 


NGC 2403 


2.1 


2 


Vi 


1.30 


1.62 


NGC 3198 


1.1 


2 


Vi 


1.75 


1.94 


NGC 0598 


1.2 


3 


Vi 


0.73 


2.02 


NGC 6503 


2.2 


3 


Vi 


1.00 


2.08 


NGC 0247 


1.1 


2 


v 2 


0.85 


3.54 


NGC 0300 


1.0 


3 


Vi 


0.62 


2.58 


347-g33 


1.2 


3 




1.66 


0.99 


UGC 2259 


1.0 


2 


v 2 


0.39 


3.9 


NGC 7339 


1.2 


2 


Vi 


1.62 


1.39 



are poorly determined in that range of distances and that the law of gravity must be assumed to be Newton's at short 
scales, allows one to discard this range in our phenomenological study. 

We have selected a sample of nine galaxies such that the luminous matter inside them can be well described by 
a thin disk with exponential density. We have applied our equation to the rotation curve of each of these galaxies 
and have found a law of gravity that can generate the observed curves without the need for dark matter (or at least, 
with a moderate quantity of dark matter distributed with the same exponential law as the luminous matter). These 
"nine laws" are statistically compatible among themselves, and point in the direction that a single, non-Newtonian, 
universal (i.e. the same for all the galaxies) law may be at work in the realm of the galaxies. 
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